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Extensive Monte Carlo simulations are performed on a two-dimensional random field Ising model. 
The purpose of the present work is to study the disorder-induced changes in the properties of 
disordered spin systems. The time evolution of the domain growth, the order parameter and the 
spin-spin correlation functions are studied in the non equilibrium regime. The dynamical evolution 
of the order parameter and the domain growth shows a power law scaling with disorder-dependent 
exponents. It is observed that, except for very small random fields, exchange interaction never wins 
over pinning interaction to establish long range order. 
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I. INTRODUCTION 

The random field Ising model (RFIM) belongs to a 
class of disordered spin models in which the disorder is 
coupled to the order parameter of the system. A lot has 
been studied on various aspects of RFIM since Imry and 
Ma [H introduced this model. It has experimental real- 
izations in diluted antiferromagnet 0]. Its Hamiltonian 
differs from that of normal Ising model by the addition 
of a local random field term which results in a drastic 
change of its behaviour in both equilibrium and non equi- 
librium situations. In one dimension {d = 1), the RFIM 
does not order at all Q; in = 2, the RFIM does not 
actually posses a order-disorder phase transition and the 
lower critical dimension di =2 [J, 0, HI and for d = 3, 
it does have a transition and it is ferromagnetic below 
its critical temperature Q. Recently Q gave numerical 
evidence that the 2d non equilibrium zero-temperature 
RFIM exhibits a critical behaviour. The Hamiltonian 
for such a system is, in general, given by 



H 



E 



E^ 



(1) 



where J is the coupling constant, conventionally set to 
unity in the present work, rji is the quenched random field 
and Hext is the external magnetic field. In the present 
work the external magnetic field H^xt is set to zero. 

The local static random fields give rise to many local 
minima of the free energy and the complexity of the ran- 
dom field free energy also gives rise to long relaxation 
times as the system lingers in a succession of local min- 
ima on its way to the lowest energy state. For the zero 
field Ising model in two or more dimensions below the 
critical temperature Tc, there forms domains of moder- 
ate size in the early time regime, in which all the spins are 
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either up or down and as time progresses the smallest of 
these domains shrink and vanish, closely followed by the 
next smallest, until eventually most of the spins on the 
lattice are pointing in the same direction. The reason 
for this behaviour is that the domains of spins possess 
a surface energy - having a domain wall cost an energy 
which increases with the length of the wall, because the 
spins on either side of the wall are pointing in opposite 
directions. The system can therefore lower its energy by 
fiipping the spins around the edge of a domain to make it 
smaller. Thus the domains "evaporate" , leaving us most 
of the spins either up or down. 

However, the story is different for RFIM. In the RFIM, 
domains still form in the ferromagnetic regime, and there 
is still a surface energy associated with the domain walls 
but it is no longer always possible to shrink the domains 
to reduce this energy. The random field acting on each 
spin in the RFIM means that it has a preferred direction. 
Furthermore, at some sites there will be a very large local 
field r]i pointing in one direction, say up direction, which 
means that the corresponding spin will really want to 
point up and it will cost the system a great deal of energy 
if it is pointing down. The acceptance ratio for fiipping 
this spin contains a factor ex'p{—2(3\r]i\), which is a very 
small number if \r]i\ is large. It is said that the domain 
wall is pinned by the local field, i.e., it is prevented from 
moving by an energy barrier produced by the large ran- 
dom field. If the domains eventually stop growing, the 
system will be in a disordered phase (although the do- 
main size may be very large). This describes the d = 2 
RFIM. 

All these features can be visualized from some snap- 
shots presented in Fig. [T] In the early time regime, 
domains begin to form and grow in size (Fig. [1] (a) - 
(d)). This is seen in pure system too. As time flows, 
the domains stop growing and the domain wall is pinned 
(Fig. [1] (e) - (h)) even at a low temperature. The intro- 
duction of a local static random field term changes the 
behaviour of the system completely in later time regime 
and renders the system in a disordered phase even at a 
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FIG. 1. (Color online) Time evolution of domains for a 128 x 128 system at T = 0.50. The strength of disorder is rjo = 1.0. 
The yellow (grey) regions indicate domains of up spins. 



3 



very low temperature. 

These disordered spin systems have been an active area 
of research for quite some time now The pur- 

pose of the present work is to study the disorder induced 
changes in the properties of these systems, with empha- 
sis on the dynamical evolution of the domain growth, the 
order parameter and the spin-spin correlation functions 
in the non equilibrium region. 

In most of the studies on RFIM, the domain size (or the 
cluster size) was determined in terms of the fluctuations 
of the magnetization [l5l - [l9j . The fluctuations in magne- 
tization is only a measure of domain size, not the actual 
domain size. In this work, we have measured the actual 
domain size by Hoshen-Kopelman algorithm [23|. The 
method of determining cluster size by Hoshen-Kopelman 
algorithm allows us to make a more accurate study of 
its growth with time. Moreover, there are many ways in 
which the random flelds could be chosen. In most of the 
studies, they were chosen to have a Gaussian distribution 
with some finite width a, or to have values randomly ±h 
where h is a constant. Although, the interesting proper- 
ties of RFIM are believed to be independent of the exact 
choice of the distribution which is a consequence of the 
phenomenon of universality [2l| . In the present work, 
the random fields are chosen from a uniform distribution 
of varying strengths. 

The rest of the paper is arranged as follows. Section 
im discusses the computational details and gives the def- 
inition of the thermodynamic quantities of our interest. 
Section IIIII presents the results in detail. Finally in Sec- 
tion |lVl we summarize our results. 

II. THE MODEL AND ITS SIMULATION 

As we have mentioned that the external magnetic field 
Hext has been set to zero in the present work, the Hamil- 
tonian reduces to 

H = -J^^SiSj + ^^i}iSi (2) 

In the present study, rji is drawn from a uniform distri- 
bution of varying strengths. The first term of the Hamil- 
tonian ^ is the usual exchange interaction term while 
the second term represents the interaction of the random 
field with the spin. We call this term the pinning inter- 
action term since this term is responsible for domain wall 
pinning. We have performed the simulations at a finite 
temperature T = 0.50, well below the critical tempera- 
ture of 2d zero field Ising model and on a square lattice 
of size 256 x 256 which can be considered as a reasonably 
large lattice size for the purpose of our investigation. We 
have chosen Metropolis algorithm |22| to simulate the 
system. Periodic boundary conditions (PBC) are used in 
the simulations. 

The sizes of the domains (or clusters) are determined 
by Hoshen-Kopelman (HK) algorithm [20]. The general 
idea of the HK algorithm is that we scan through the 



lattice looking for up spins (or down spins). To each up 
spin (or down spin) we assign a label corresponding to the 
cluster to which the up spin (or the down spin) belongs. 
If the up spin (or the down spin) has zero neighbours of 
same sign, then we assign to it a cluster label we have not 
yet used (it's a new cluster). If the up spin (or the down 
spin) has one neighbour of same sign, then we assign to 
the current spin the same label as the previous spin (they 
are part of the same cluster) . If the up spin (or the down 
spin) has more than one neighbour of same sign, then 
we choose the lowest-numbered cluster label of the up 
spins (or the down spins) to use the label for the current 
spin. Furthermore, if these neighbouring spins have dif- 
ferent labels, we must make a note that these different 
labels correspond to the same cluster. The HK algorithm 
is a very efficient cluster identification method for two- 
dimensional systems. The domain size corresponds to the 
number of spins enclosed by the boundary of a domain. 

The time dependent ensemble average of magnetiza- 
tion, i.e., the order parameter has been defined as 

Af(,7o,t) = (-^|^s, |)t (3) 

i 

where L is the linear size of the system and Si is the spin 
at the site i. The angular bracket {■ ■ ■)t indicates the 
ensemble average at time t. 

The time dependent ensemble average of the spin-spin 
correlation function is defined as 

ipssivoJjt) = {{siSi+i)i)t (4) 

where I is the distance of separation between the spins. 
The angular bracket (••■); indicates the ensemble average 
over the distance of separation between the spins while 
(• • • )t indicates the same over time t. 

The thermodynamic quantities of our interest are av- 
eraged over 50 independent simulations to improve the 
accuracy and the quality of the results. 

III. RESULTS AND DISCUSSIONS 

The dynamic evolution of the order parameter defined 
in Eq. ([3]) for different disorder strengths is shown in 
FIG. m It is seen from FIG. [2] that for weak disorder 
ivo ^ 0.3) the system reaches in a steady state after a 
certain time (here steady state means the fluctuations 
in order parameter is very small with time). For weak 
disorder, the system behaves more or less like a pure sys- 
tem. For larger disorder, the system takes longer time 
to reach in a steady state, i.e., the system relaxes, if at 
all, very slowly. The nature of the evolution of the or- 
der parameter strongly depends on the strength of the 
disorder. The early time behaviour of the dynamic evo- 
lution of the order parameter can be characterized by the 
following power law behaviour 

M(77o,i)~i^^"°^ (5) 
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exchange interaction defined as 



(6) 



as well as the time evolution of the time-dependent en- 
semble average of the normalized pinning interaction de- 
fined as 



^{Vo,t) = {—{sir],))t 



(7) 



where the angular bracket {■ ■ ■)t indicates the ensemble 
average at time t. FIG. 2] shows the dynamic evolution 
of the x{Vo,t) and the ^l{T]o,t) respectively and it reveals 



FIG. 2. (Color online) Plot of the time evolution of order 
parameter. The dashed lines represent the best linear fits 
according to the scaling law defined in Eq. ((Sj 



where /3{r]o) is a disordcr-strcngth-dcpendcnt exponent 
corresponding to the growth of the order parameter. The 
dependence of the power law exponent P{rio) on the dis- 
order strength rjQ is shown in FIG.[31 The exponent /3{rio) 
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FIG. 3. (Color online) Variation of the exponent /3(??o) with 
■qo. The solid line indicates the best exponential fit to the 
data points. 

falls off exponentially with the strength of the disorder 770 
as ~ exp(— 1/770) with ly = 1.3527 ± 0.0730. This is quite 
obvious. Two types of interactions, namely, the exchange 
interaction and the pining interaction are present in the 
system. For larger disorder, spin flips are not favoured 
due to the presence of the pinning interaction term and 
consequently the domain walls get pinned, leaving the 
system in a disordered phase. As a result, the exponent 
13{t]o) falls off sharply with 770. 

What happens is that there is always a competition 
between the exchange interaction and the pining inter- 
action. Therefore it is interesting to observe the time 
evolution of the time-dependent ensemble average of the 
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FIG. 4. (Color online) The time evolution of the exchange 
interaction xivojt) and the pinning interaction fi(r;o, t) as de- 
fined in Eq. ([6]) and Eq. ((7| respectively. 

a striking feature. It is seen that the exchange interac- 
tion xivojt) remains almost same with time except at 
some initial time steps whereas the pinning interaction 
ri(77o,t) decays more rapidly except for very large dis- 
order strength. The system evolves with time in such 
a way that the pining interaction gets minimized. We 
would like to point out that this behaviour is character- 
istic to the RFIM and was not observed earlier to the 
best of our knowledge. For weak disorder, the system 
achieves the steady state with the flow of time resulting 
in the saturation of the pining interaction and the ex- 
change interaction as well. For disorder strengths lying 
in the intermediate range, the decay of pining interaction 
strongly depends on the strength of the disorder. Now a 
natural question that arises is that how rapidly does the 
pining interaction 0(770, ^) decay or how slowly does the 
exchange interaction xivoit) increase with the strength 
of the random field 770. According to the nature of the 
dynamic behaviour seen in FIG.|4l the power law depen- 
dence of the xivoit) and the Q{riQ,t) with time has been 
proposed as 

xivo,t) ^ (8) 
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FIG. 5. (Color online) Variation of 0^(770) and an{rio) with 770- 



where a^^ijo) and an(?7o) are two characteristic expo- 
nents which determine the nature of growth of the ex- 
change interaction and the pining interaction. The vari- 
ation of these exponents with r/o is shown in FIG. [S] 
The values of the exponent a~t,(r]o) are very small and 
its variation with the strength of the disorder 770 is neg- 
ligible which indicate that 0(^{riQ) is almost independent 
of 770. On the other hand, the variation of the exponent 
ctnivo) with tjq is noticeable and 00(770) decreases rapidly 
with riQ and tends to saturate for very large values of rjo ■ 
These results confirm our earlier observation on the time 
evolution of the exchange interaction and the pinning in- 
teraction. 

Next we focus our attention on the study of spin-spin 
correlation functions defined by Eq. The spin-spin 
correlation functions tpssivoi^) for different strengths of 
disorder at different times t are shown in FIG. [51 It is 
seen from FIG.IHlthat the system is being correlated with 
time. For weak disorder, the system possesses long-range 
order after long time evolution. As time progresses the 
pinning interaction gets minimized and then saturates 
for weak disorder while the exchange interaction remains 
almost unchanged. Thus short-range order at the early 
time regime and long range order at later time regime 
are being established in the system. This is why, we ob- 
tained a large value of the exponents /3(?7o) and 0:0(770) 
for weak disorder strengths. For strengths of disorder ly- 
ing in the intermediate range, no long range order is seen 
to develop and short range order is being developed with 
time in the system. Exchange interaction never wins over 
pinning interaction to establish long range order in the 
system. Due to this short-range order, there forms small 
domains in the system initially and as time elapses, these 
small domains "evaporate" to form a large domain. How- 
ever, the presence of random fields prevents the system 
from growing into a single domain even at a low tem- 
perature. For very large disorder (> 2.6), neither short 
range order nor long range order prevails in the system 



and the system remains in a complete disordered phase. 

Now we present the results of the domain growth with 
time. Domain growth in quenched non equilibrium sys- 
tems is a widely studied topic I2J . Monte Carlo simu- 
lations have been carried out |24l427| to obtain numeri- 
cally the time evolution of these domains, and to compare 
the results with theory j2§ - [3]| . We concentrate on the 
growth of the largest domain and it is plotted in FIG. [T] 
The growth of the largest domain can be characterized 
by the following power law behaviour : 
For 7/0 < 0.3 



A(?/o,t) = 




t < txl 


Mvo,t) = 










t > tx2 



For 770 > 0.3 

A(7;o,i) = t^^''°^ (10) 

where /i(?7o) is a disordcr-strcngth-depcndent exponent 
corresponding to the growth of the largest domain at 
the early time regime. The variation of the exponent 
fj'ivo) with 770 is shown in FIG. |8l The exponent iJ-{rio) 
falls off exponentially with the strength of disorder 770 
as - exp(-777o) with 7 = 1.9276 ± 0.1409. It is to 
be noted that although the order parameter exponent 
and the growth exponent falls off exponentially with the 
strength of disorder, the values of the exponents are dif- 
ferent and the fall of A<(?7o) is faster than that of /3(?7o)- 



IV. SUMMARY AND CONCLUSION 

We conclude the paper with a summary of our re- 
sults. This paper has attempted to consider some as- 
pects of the non equilibrium behaviour of the 2d random- 
field Ising model numerically at a low temperature. As 
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FIG. 6. (Color online) Plot of ipssivoj) against I for different disorder strengths at different times. 




FIG. 7. (Color online) Plot of the growth of the largest do- 
main of the system with time. The dotted lines represent the 
best linear fit to the data points. 
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FIG. 8. (Color online) Plot of the exponents /i(7?o) with dis- 
order strength r/o 



seen in the preceding Sections, the RFIM exhibits a va- 
riety of behaviour depending on the strength of the ran- 
dom field. The dynamic evolution of order parameter 



shows a power-law growth at the early time regime with 
a disorder-dependent exponent, we find from the study 
of spin-spin correlation functions that the RFIM shows 
short-range order, long-range order and no order at all. 
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each of which occurs in a restricted range of random 
field strength. This observation suggests the possibil- 
ity of existence of a critical field strength up to which 
the domain size can grow until it reaches a maximum. 
However, a more detailed study is required to make a 
definite statement in this regard. We mention in pass- 
ing that [3^ has shown the existence of critical Gaussian 
random field above which the geometrical clusters have 
a finite extent and below which they are percolating. It 
was shown that the system evolves with time in such a 
way that the pinning interaction gets minimized whereas 
the exchange interaction remains almost unchanged. Ex- 
cept for very weak disorder, exchange interaction never 
wins over pinning interaction to establish long range or- 
der. This is because a small random field r]i can be 
looked upon as another "spin" in the environment of Si, 



rather than a pinning force. The domain growth in this 
quenched non equilibrium system reveals that its relax- 
ation breaks into several distinct time regimes for small 
random fields. Above a threshold, the dynamic evolution 
of domain exhibits a power-law growth with a disorder- 
dependent exponent. Recently, Corberi et. al. (33l | also 
found that the domain growth shows a power law scal- 
ing with a disorder dependent exponent in preasymptotic 
regime. 
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